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Abstract 


A series of analyses have been developed which permit the calculation of the performance of 
common inlet designs. The methods presented are useful for determining the inlet weight 
flows, total pressure recovery, and aerodynamic drag coefficients for given inlet geometric 
designs. Limited geometric input data is required to use this inlet performance prediction 
methodology. The analyses presented here may also be used to perform inlet preliminary 
design studies. The calculated inlet performance parameters may be used in subsequent 
engine cycle analyses or installed engine performance calculations for existing uninstalled 
engine data. 


Introduction 

Propulsion installations can have a significant effect on the overall efficiency of airbreathing 
engine systems, particularly for supersonic and hypersonic flight vehicles. To assess the 
impact of an inlet design on the net thrust and specific fuel consumption for a given engine 
design, either the inlet performance characteristics must be known in advance, or they must 
be calculated from a simple geometric design, or in the worst case the inlet system must be 
designed from scratch and then analyzed to determine performance. This report describes a 
series of analyses which have been developed into a performance prediction methodology for 
engine inlet systems. The methodology can be used to predict performance for a given inlet 
geometric design. Additionally, the methodology can be employed to perform preliminary 
inlet system design, and subsequent performance analyses. 

Inlet performance is typically comprised by determining three quantities: delivered engine 
airflow, W 2 , total pressure recovery, Pn/P m , and aerodynamic drag coefficient, C D . It is 
also very important to be able to characterize inlet performance over the entire vehicle flight 
and engine operation range, not just at the inlet design point. The methodology presented 
covers the calculation procedures used to determine inlet performance, both on and off- 
design, for the broad classification of inlet geometries shown in Figure 1. 
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The geometric input used for the analysis modeling is simple and flexible. This permits 
rapid performance calculations and quick turn-around times for inlet design assessments. 
The analyses are capable of modeling three broad inlet design classifications: pitot, 
axisymmetric, and two-dimensional. 


Method of Analysis 

Figure 2 shows the basic modeling elements used to develop the inlet performance analysis 
methodology. The action of airflow ingestion through the inlet is broken up into a series of 
distinct processes. Changes in flow properties from the free stream flow station, 0, to the 
inlet local flow station, L, are modeled as vehicle effects. Flow changes through shock 
waves ahead of the cowl lip station, 1, are modeled as external compression. Flow changes 
within the cowl lip to the inlet throat station, TH, are modeled as internal compression. 

Flow changes downstream of the throat to the engine face station, 2, are modeled as subsonic 
diffusion. 

Aerodynamic drags modeled include spillage, bleed, and bypass. Spillage drag is the sum of 
the momentum change incurred by air being diverted around the inlet lip, additive drag, and 
cowl lip suction, if present. Bleed drag results from the momentum change in air which is 
dumped overboard as required by inlet stability considerations and boundary layer control. 
Bypass drag results from the momentum change in air which is dumped overboard for 
inlet/engine weight flow matching requirements. Additional calculations for cowl lip and 
wave drag are also included. 

The relative amounts of airflow ingested into the inlet, lost to bleed, bypass, or spillage are 
shown in Figure 2 as the free stream tube areas, A. These areas are usually presented in 
analyses as a ratio with respect to the forward projected cowl lip area, A c . 

Figure 3 shows the different modes of operation which are possible for mixed-compression 
inlets in supersonic flight. Of particular importance is the location of the normal shock wave 
since this will dramatically effect the inlet airflow capture characteristics. The top diagram 
in Figure 3 shows the inlet operating with the normal shock wave outside of the cowl lip. In 
this mode of operation the inlet can deliver less (or within limits more) air to the engine by 
spilling air around the lip as subsonic flow behind the normal shock wave. Thus the engine 
demand can influence the inlet operation and the location of the normal shock wave. This 
operation is called sub-critical and the inlet is unstarted. External compression inlets always 
operate sub-critical. 

As the engine demands more airflow, the normal shock wave is drawn up to the cowl lip. 
When the normal shock wave just reaches the cowl lip, the inlet is ingesting the ma vinmni 
airflow possible. The center diagram in Figure 3 shows this operation, called critical, but 
the inlet is still unstarted since the throat Mach number is subsonic. When the normal shock 
wave is swallowed and located downstream of the throat the operation is called super-critical 
and the inlet is now started since the throat Mach number is supersonic. The airflow 
captured by the inlet lip cannot be increased or decreased by the engine operation, but is 
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fixed by the external shock wave structure as shown in the bottom diagram of Figure 3. 
Inl et/engine airflow matching can only be accomplished in this mode using a bypass system. 


Engine Weight Flows 

The primary function of an inlet is to deliver the proper amount of airflow to the engine. 

The amount of airflow delivered to the engine depends on many factors. The usual 
requirements for inlet design specify the desired altitude corrected weight flow delivered to 
the engine face as a function of flight Mach number. Equation 1 shows the relation between 
engine corrected weight flow and inlet performance and design variables. The leading term 
in Equation 1 is the inlet capture area, A c . The larger the inlet the greater the engine weight 
flow. The second term is the free stream- tube area ratio and it is a strong function of the 
inlet design and mode of operation. Stream-tube area ratios will be discussed in a later 
section. 

The inlet total pressure recovery is also an important factor in Equation 1 , however, the 
corrected airflow is inversely proportional to recovery. A higher recovery will result in a 
lower specific corrected airflow at the engine face, and hence will necessitate a larger inlet, 
A c . And this in turn will result in a propulsion system capturing more absolute airflow, 
resulting in greater thrust. A lower recovery results in a smaller inlet, less absolute airflow, 
and lower thrust. 

Equations 2 through 8 show how the absolute engine weight flow, W 2 , is calculated from the 
corrected weight flow. Equation 9 indicates that the free stream static pressure and 
temperature are known from standard atmosphere tables or curve fits. Equations 7 and 8 are 
the isentropic flow relations between total and static quantities as a function of Mach number. 
Equation 6 is a statement of the first law of thermodynamics, and is only valid if the inlet 
does not transfer heat or shaft work to or from the airflow. Equation 10 may be used to 
determine the actual weight flow directly from free stream static properties. Equations 1 
through 10, as written, imply the use of English units. 


Real Gas Effects 


Implicit in Equations 1 through 8 is the ideal gas assumption. This is usually valid for free 
stream Mach numbers below two. At higher flight speeds real gas effects need to be 
accounted for. Equations 11 through 15 are used for a calorically imperfect gas model. 
Primed values correspond to the real gas property. Equation 11 is used to calculate a real 
gas ratio of specific heats from the ideal gas 7 and the static temperature. For a known 
flight Mach number Equation 12 is solved by iteration to yield a real gas total temperature. 
Equation 13 is used to determine the total pressure for the real gas model. These total 
quantities are used to replace the ideal gas values calculated by Equations 7 and 8 . Equation 
1 must also be modified for stream-tube area variations using Equation 14. Additional 
information on this real gas model can be found in reference 1 . 
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Inlet Mass Flow Ratios 


Figure 2 shows a rather standard airflow accounting system in terms of idealized free stream- 
tube areas. If these stream-tube areas are normalized by the inlet capture area, a series of 
relations can be developed. Equations 16 through 18 show the stream-tube area build-up. 
These ratios are also called mass flow ratios, since the mass flow is equal to density times 
velocity times area. The density and velocity terms drop out in the ratio format. If there are 
no vehicle effects on the airflow ahead of the inlet, the right hand side of Equation 18 is 
equal to one. The airflow captured by the cowl lip and ingested into the inlet is represented 
by the mass flow ratio A 0I /A C , while the airflow passed through the inlet throat is represented 
by the mass flow ratio Ag/A c . 


Vehicle Effects 


The effects of a vehicle flow field ahead of an inlet can be simply described as chang es to 
the total pressure, Mach number, and stream-tube area between stations 0 and L. Equations 
19 through 25 show the effects of Mach number and total pressure changes on the stream 
tube area. The analysis extends from the principle of conservation of mass in Equation 19. 
Equation 25 provides a simple expression for determining the right hand side of Equation 18 
if the total pressure ratio and Mach number ratio are known from stations 0 to L. 

All of the subsequent analyses are performed in the inlet local reference frame, as if there 
were no vehicle effects present and the inlet was simply in a free stream of different Mach 
number and total pressure. However, the overall inlet performance must be represented in 
the free stream reference, and thus all results from the inlet local reference must be adjusted. 
Equations 26 through 31 show how a drag coefficient calculated in the inlet local reference is 
adjusted to represent the same force described as a drag coefficient in the free stream 
reference frame. 


Vehicle Forebodv Model 

Figure 4 shows a simple vehicle forebody model employed in the methodology. This model 
can be used to represent vehicle underbody precompression surfaces, upperbody expansion 
surfaces, aircraft wings, or slender fuselages. The stream tube area shown in Figure 4 can 
be seen to decrease as the streamline crosses subsequent shock waves from the free stream to 
the inlet local stations. Equations 32 and 33 show how the ratios of total pressure and Mach 
number, from free stream to inlet local, are determined from the changes across each 
individual flow deflection region. Positive angles, a, are modeled as oblique shock wave 
compression regions. Negative angles are modeled as discrete Prandtl-Meyer expansion 
regions. Conic shock waves are also an included option, in addition to the default planar 
shock wave calculations. Since shock wave calculations are a central part of the inlet 
performance methodology, a description of these types of calculations follows. 
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Normal Shock Wave Relations 


A shock wave is a very thin layer interaction between two distinct compressible flow regions. 
The simplest shock wave type is the normal shock wave, shown in the top diagram of Figure 
5. Supersonic flow is shocked down to subsonic flow across a normal shock wave. All flow 
properties are determined by the upstream conditions. Equations 34-38 show the standard 
normal shock wave flow relations for Mach number, pressure, temperature, and density. 

Note that all of these relations are only a function of upstream Mach number, and are thus 
easy to apply. 


Oblique Shock Wave Relations 

The center diagram in Figure 5 shows the elements of an oblique shock wave. A planar 
oblique shock wave is produced by a downstream boundary turning an incoming supersonic 
flow through an angle 0. As a result, a shock wave forms, inclined to the incoming flow 
direction at an angle 0. The components of the flow perpendicular to the oblique shock 
wave are described by the normal shock wave relations. The velocity components parallel to 
the shock wave are unchanged. Equation 39 gives a relation between the flow turning and 
shock wave angles. Typically the flow deflection angle is known, and the shock wave angle 
must be found. Although most references suggest solving Equation 39 iteratively, this is not 
necessary. 

By some algebraic manipulation, Equation 39 can be rewritten in the form of a 6th order 
polynomial in terms of the sine of the shock wave angle, shown in Equation 40. The 
coefficient terms of the resulting polynomial are given in Equations 41 through 43. Since 
Equation 40 only has even power terms, a generalized solution for 3rd order polynomials can 
be employed. This gives a relation for the square of the sine of the shock wave angle as a 
function of the flow deflection angle and the upstream Mach number. Equations 44 and 45 
show this direct solution. Once the shock wave angle is known, the normal components of 
the Mach number are found in Equations 46 and 47. The changes in flow properties are then 
calculated by application of the normal shock wave relations. 


Conical Supersonic Flow Relations 

The calc ulati ons involved in determining supersonic flow in conical shock fields are a bit 
more complex. The elements in the conical shock wave problem are shown in the bottom 
diagram of Figure 5. Equation 48 gives the reduced differential equation describing the flow 
field between the conical shock wave and the cone surface. The dependent variable in 
Equation 48, V/, is the radial component of a non-dimensionalized velocity in the conic flow 
field. Equations 49 through 52 give definitions of the non-dimensional velocity components 
and their relation to the polar angle <f> and flow direction angle 0. Equations 53 and 54 are 
the two required boundary conditions of tangent flow to the cone surface and the correct flow 
fuming angle behind an oblique shock wave. 
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Equation 48 can be solved numerically by a scheme commonly known as the Taylor-Maccoll 
solution. A conic shock wave angle is first guessed, and with Equation 54 provides a value 
for the shock wave boundary condition. Equation 48 is then solved at small increments of <j) 
using standard Runge-Kutta integration schemes. The solution is then marched by <£ through 
the flow field to the cone angle 0 C . If the tangent flow boundary condition. Equation 53 is 
satisfied, then the initial guess on the conic shock wave angle is correct. Otherwise another 
guess on the angle /3 is chosen, and the process is repeated, iterating to a correct solution. 
Once solved, the flow properties across the conic shock wave are determined from the 
oblique shock wave relations. Also, the flow velocity variations from the conic shock wave 
to the cone surface are known from the solution of Equation 48. Other flow properties can 
be then determined from the isentropic flow relations given below. Reference 2 is a good 
starting point for further information on calculating conical shock waves. 


Isentropic Flow Relations 

Equations 55 through 59 are a series of often used isentropic flow relations found throughout 
the methodology, and are given here for convenience. Equations 55 through 57 calculate the 
static pressure, temperature, and density as functions of Mach number only. For example, 
the static pressure field behind the conic shock wave is determined by Equation 56, since the 
Mach number field is known from the solution to Equation 48 and the total pressure is a 
constant, whose value is determined from the oblique shock wave relations. Equation 58 
describes the required stream-tube flow area as a function of Mach number, where A. is the 
flow area at the sonic condition. 

Equation 59 is the Prandtl-Meyer function and it is used to determine the change in Mach 
number as a supersonic flow isentropically expands through a turning angle. Typically an 
initial Mach number is known as well as the turning angle. Equation 59 determines the 
initial Prandtl-Meyer function value explicitly. By adding the expansion angle (in radians) a 
new Prandtl-Meyer function value is calculated, from which Equation 59 must be solved 
iteratively to yield a new value of the Mach number downstream of the expansion. This 
technique is used as part of the vehicle forebody model for supersonic flow expansions. 


Total Pressure Recovery 


The total pressure recovery for the entire inlet is calculated as the product of a series of total 
pressure ratios across elements of the inlet system. Equation 60 shows this relation, where 
the terms on the right hand side are the total pressure ratios from: free stream to inlet local, 
inlet local to inlet lip, inlet lip to throat, and throat to engine face. Each of these terms is 
calculated in the subsequent modeling elements, with the exception of the free stream to inlet 
local term, which is calculated in the vehicle forebody model previously discussed. 
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External Compression 


The changes in flow properties from the inlet local station to the inlet lip are determined by 
models of the external compression processes for a given inlet design. Figure 6 shows the 
elements of the external compression models used in the methodology. Each basic inlet type 
must be modeled separately, since the external flow is highly dependent on the inlet 
geometry. 

The top diagram in Figure 6 shows the elements of the external compression model for pitot 
inlets. The total pressure ratio from inlet local to inlet lip is given in Equation 61 and the 
total pressure loss is only generated by a normal shock wave at the inlet local Mach number. 
If the inlet local Mach number is subsonic, then the total pressure ratio is one. There is an 
incurred drag penalty for air which is spilled around the cowl lip called additive drag. 
Reference 3 gives a procedure for calculating this drag coefficient and Equations 62 through 
64 summarize the analysis. The last term in Equation 62 is the mass flow ratio ingested by 
the inlet lip, A u /A c , and this number is determined by the engine airflow requirements. 
Equations 63 and 64 result from conservation of mass and the isentropic flow functions. 

The center diagram in Figure 6 shows the elements of the external compression model for 
axisymmetric inlets. This type of inlet is capable of operating either super-critical or sub- 
critical, and the model must distinguish the difference. Equations 65-67 pertain to the super- 
critical operation mode. The total pressure ratio is produced entirely by the inlet conic shock 
wave. The additive drag coefficient can be determined either by Equation 66 or 67 . 

However, since the conic flow field is known for supersonic operation. Equation 67 is 
employed using numerical integration techniques. The integration path corresponds to the 
streamline intersecting the cowl lip. For subsonic flows Equation 66 must be used and the 
total pressure ratio is one. 

For sub-critical operation, a normal shock wave exists outside of the inlet cowl lip. This 
results in a greater pressure loss, higher additive drag coefficient, and lower mass flow ratio. 
Equations 68 through 74 show these calculations for sub-critical inlet operation. The position 
of the normal shock wave outside of the cowl lip is approximated as standoff distance which 
is proportional to the inlet capture mass flow ratio relative to critical operation, as indicated 
by Equation 73. The proportionality factor, K, is a function of Mach number (indicated by 
Equation 74) and this function was determined from curve fits to data found in reference 4. 
The functional form of the shock wave standoff factor, K, is shown graphically in Figure 7. 

The bottom diagram in Figure 6 shows the elements of the external compression model for 
multi-ramp two-dimensional inlets. These inlet types can also operate both sub-critical and 
super-critical. Equations 75-78 show the calculations for super-critical operation. The total 
pressure ratio in Equation 75 is the product of all the external oblique shock wave total 
pressure ratios. For sub-critical operation, Equations 79-82 show the calculations used in the 
model. A normal shock wave can exist outside of the cowl lip and the relations computing 
the total pressure ratio and additive drag need to account for the position of the normal shock 
wave, and on which ramp it is located. 
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Internal Compression 


Figure 8 shows the elements of the internal compression model. For started inlet operation, 
an oblique shock wave train is used to model the losses in the internal portion of the inlet 
from the cowl lip to the throat. The net turning angle is the sum of the last external surface 
angle and the internal cowl lip angle. The flow properties across each shock wave reflection 
are determined from the oblique shock wave relations previously discussed. Equation 83 
shows the relation between the flow properties in the model and the geometric throat area 
constraint. Equation 83 is again a statement of conservation of mass for compressible flows. 

The reflecting oblique shock wave model, Equations 85 through 87, is primarily used to 
determine the total pressure loss in the internal compression region. The model may also be 
used to determine a throat Mach number, M m , for a given throat area ratio, A m /A c , by 
iterative solution of Equation 83. Often the throat Mach number is specified instead and the 
throat area ratio is then determined directly by Equation 83. If both the throat Mach number 
and throat area ratio are specified, then the inlet capture mass flow ratio, AJA C , must then 
be determined from these constraints. 


Subsonic Diffusion 

Figure 9 shows the elements of the subsonic diffusion model. Depending on the inlet 
operation mode, a terminal normal shock wave may or may not exist downstream of the inlet 
throat within the subsonic diffuser. Equations 88 through 94 show the calculation procedure 
for operation with subsonic flow at the inlet throat. The model used here closely follows that 
given in reference 5. For inlet operation with a subsonic throat, the throat Mach number and 
area are usually specified, consistent with the desired engine weight flow delivered. For 
inlet operation with a supersonic throat, or started operation. Equations 95 through 97 are 
used. The strength of the terminal normal shock wave can be used to provide inlet/engine 
corrected mass flow matching in some instances. Curve fits are used for the loss factor 
functions in Equations 90 and 94 corresponding to divergence and throat Mach number loss 
mechanisms. Figures 10 and 11 show the functional forms for the divergence and throat 
Mach number loss mechanisms graphically. The friction factor given in Equation 93 is a 
nominal value, and may be changed if desired. 


Bleed Drag 

Bleed drag seems to be a necessary evil required for supersonic inlet designs. Since inlets 
produce large positive pressure gradients, some severe in shock wave interactions, the 
boundary layers are prone to separation. To alleviate this problem, portions of the boundary 
layer are removed through wall suction, and then dumped overboard. If done correctly, this 
usually results in improved inlet recoveries, however, a momentum drag is incurred. 
Equations 98 through 109 show the procedure used to calculate the bleed drag coefficient. 
These relations follow the procedures outlined in reference 6. Total pressure losses up to the 
bleed system plenum are modeled, as well as the effective bleed nozzle exit pressure and 
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flow area. Non-axial nozzle exit flow losses are also included. 

A number of inputs must be specified for the design and operation of the bleed system. In a 
high speed inlet the bleed system is typically comprised of a series of discrete bleed regions, 
each having its own type of wall perforation, plenum, and exhaust nozzle. The total bleed 
drag is thus the sum of the individual bleed system elements. Equation 98 is used to 
describe the bleed drag for a discrete bleed element. The bleed mass flow ratio, A LBLD /A C , 
nozzle exhaust flow angle, 6 X , and nozzle exhaust velocity coefficient, tj v , must all be 
specified. Additionally, the bleed plenum recovery, PtrJPtl , must also be specified. To 
choose these values extensive experience in the design and operation of bleed systems is 
usually required. To alleviate this requirement default values have been implemented in the 
methodology. 

Figure 12 shows typical bleed system operating characteristics which are incorporated as user 
selectable defaults for inputs to the bleed drag model. The bleed plenum recovery is shown 
as a function of inlet local Mach number for a variety of bleed system design wall 
perforations. The total bleed mass flow ratio required for typical inlet operation is also 
shown in Figure 12 over the same Mach number range. The data which comprises the basis 
for Figure 12 is taken from reference 7. 

There are two additional empirical relations embedded within the bleed drag model. 

Equation 99 shows the functional dependence for the oblique exit nozzle drag factor, C re , 
and Figure 13 shows this functional relationship graphically. The relationship for the 
effective nozzle discharge pressure, Equation 103, is shown graphically in Figure 14. The 
bleed exhaust nozzle area ratio, A x /A w , is the final input required for the bleed drag model. 
The nozzle area ratio should be chosen depending on the bleed exhaust nozzle pressure ratio. 
The operating pressure ratio for the bleed exhaust nozzle is given in Equation 102. Based on 
this value, Figure 15 can be used to pick the appropriate nozzle area ratio for the bleed 
element. Other area ratio choices will result in over or under expansion losses which will 
further increase the resulting bleed drag. Since bleed plenum recoveries are typically low 
this usually results in the use of convergent nozzles for bleed systems. Bypass systems can 
have much higher recoveries, and thus may be able to utilize convergent-divergent nozzle 
designs. 


Bypass Drag 

Bypass flows are used to dump air overboard in the subsonic diffuser ahead of the engine 
face, and are typically employed for inlet/engine flow matching. The resulting drag 
coefficient is calculated in a mann er analogous to that used for the bleed system. Equations 
110 through 115 show the modifications made to the bleed drag relations required to model 
bypass flow. The required inputs to the bypass drag model parallel those necessary for the 
bleed drag model. As in the bleed system model, the bypass system can be comprised from 
a series of distinct bypass elements. Each element can be defined with different design and 
performance characteristics. The total bypass drag thus being the sum of the drags of all the 
distinct elements. The bypass plenum recovery is typically a function of the amount of 
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bypass flow dumped overboard. Figure 16 shows the methodology default for the relation of 
the bypass recovery, P^p/Pn, as the bypass mass flow ratio, A^yjJAq, varies. The data on 
which Figure 16 is based can be found from reference 7. 


Cowl Lip Suction 

As a result of sub-critical airflow spillage around the inlet cowl lip, the static pressure over 
the cowl leading edges is decreased, thus reducing the effective cowl pressure drag. This 
effect, known as cowl lip suction, can be viewed as a correction to the additive drag 
calculation as presented previously in the external compression model. The net combination 
of the additive drag and cowl lip suction is the total inlet spillage drag. Equation 116 shows 
the definition of the cowl lip suction coefficient. This model for the cowl lip suction 
coefficient is based entirely on empirical relations which can be found in reference 6. 

Equations 117 through 124 detail the empirical terms used in Equation 116. The functional 
form of Equation 117, the first cowl lip suction factor, K a , is shown graphically in Figure 
17. The effective cowl lip angle correction factor, a, defined in Equation 118 is shown 
graphically in Figure 18. The procedure for computing the effective cowl lip angle is given 
in Equations 119 through 121. Equation 119 is an approximation for the effective cowl lip 
angle, in degrees, determined from the integral parameter, G, which is defined by Equation 

120. This integral parameter evaluates the cowl surface curvature from the cowl lip leading 
edge to the maximum of the cowl forward projected area location. In Equations 120 and 

121, the cowl profile is defined by coordinates (X,Y) and the cowl lip leading edge is located 
at (X C ,Y C ). 


The second cowl lip suction factor, K 0 , defined in Equation 122 is shown graphically in 
Figure 19. The final empirical cowl lip suction factor, C D2 , is defined in Equation 123 and 
is also shown graphically in Figure 20. Once the cowl lip suction factors are determined 
from curve fits and the cowl lip suction coefficient calculated, the inlet spillage drag 
coefficient is then found by Equation 124. 


Cowl Lip and Wave Drag 

The pressure drag acting on the inlet cowl surfaces can typically be broken into two parts; 
drag due to a blunt inlet lip and wave drag due to the area growth along the remainder of the 
cowl surfaces. Equation 125 shows this drag decomposition. For sharp lip inlets, the drag 
component due to a blunt lip is necessarily zero. For non-sharp lip inlets, the blunt leading 
edge will produce a pressure drag at supersonic local Mach numbers resulting from a 
detached normal shock wave which is formed over the leading edge radius of the cowl lip. 
Equation 126 shows the computation of the lip drag coefficient based on the assumption that 
an average pressure rise produced by a normal shock wave at the inlet local Mach number 
acts over the forward projected cowl lip surface area. This average pressure rise is modeled 
as the simple arithmetic mean of the stagnation and static pressures behind a normal shock 
wave. The forward projection of the blunt cowl lip area is denoted as A x in Equation 126. 


10 


The pressure drag acting on the rest of the cowl surface area is wave drag. The wave drag 
coefficient is defined by Equation 127 for two-dimensional inlet geometries. If the cowl 
profile is comprised by a series of flat plates, the integration in Equation 127 can be replaced 
by a discrete summation as shown in Equation 128. The pressure acting on each cowl plate 
segment, P„ is calculated by the shock wave and expansion models previously described. 

The forward facing projected area of each cowl segment plate is denoted as in Equation 
128. 

The computations of the wave drag for axisymmetric cowls are given in Equations 129 
thr ough 143. The wave drag coefficient is defined as an integration of the pressure 
coefficient over the cowl surface as shown in Equation 129. Equation 129 is an equivalent 
statement to Equation 127 which defined the wave drag coefficient for two-dimensional inlet 
geometries. The computation of the pressure coefficient, C P , over an axisymmetric cowl 
geometry, however, is substantially more complex than the two-dimensional flat plate cowl 
model of Equation 128. The pressure coefficient in the axisymmetric wave drag model, 
given in Equation 130, is calculated by a first order approximation using the perturbation 
velocities determined from the solution of a linearized supersonic slender body theory. 

The axisymmetric form of the governing partial differential equation for the perturbation 
velocity potential by supersonic slender body theory is given in Equation 131. The 
generalized solution of the perturbation velocity potential, <f>, and the axial and radial 
perturbation velocities, u and v respectively, are found in reference 8 and given in Equations 
132 through 135. In Equations 133 and 134,/ '(f) is a singularity distribution along the 
centerline axis which uniquely determines the flow field on and about the slender body 
surface. A statement that the flow is tangent to the body surface on the body surface can be 
used as a boundary condition to determine the singularity distribution for that body. A more 
detailed description of the analyses which follow can be found in reference 9. 

If the axisymmetric cowl surface profile is described by the coordinate pairs ( X,R ) then the 
body surface tangent flow boundary condition can be written as Equation 136. Furthermore, 
if the cowl surface profile is discretized and the singularity distribution, / ', can be assumed 
piece-wise constant over a small interval then the discrete elements of the singularity 

distribution can be moved outside of the integration, as shown in Equation 137. The initial 
and final bounds of the piece-wise integrations are given in Equations 138 and 139 as they 
apply to Equation 137. The piece-wise integral is now readily evaluated in closed form, and 
the solution becomes Equation 140. A marching scheme can easily be developed to 
determine the value of a discrete singularity, /„', corresponding to a location (X„,R„) on the 
cowl surface in terms of a summation of all the upstream singularities, as shown in Equation 
141. Therefore, the entire singularity distribution can be determined by simply marching 
down the cowl surface using Equation 141. 

Once the discrete singularity distribution is known, the pressure coefficient can be 
determined by an analogous procedure, as shown in Equation 142, which is developed from 
Equations 130 and 133. Again, the resulting piece- wise integral in Equation 142 can be 
evaluated in closed form, yielding Equation 143, and the pressure coefficient at discrete 
points along the cowl surface is subsequently known as a function of the discretized 
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singularity distribution. The axisymmetric wave drag coefficient is then determined by 
numerical integration of Equation 130 using the values found from Equation 143. 


Lip Losses 

The inlet cowl lip can have additional effects on the inlet recovery, particularly at low 
speeds. At take-off conditions, the inlet must ingest mass by drawing a large volume of 
initially stationary air from the surroundings around the cowl lip and then into the engine 
face. For sharp lip inlets, as the airflow is drawn around the cowl lip, the flow will 
accelerate and separate as it turns, producing a subsequent fluid dynamic loss and drop in 
total pressure recovery. Equation 144 shows the total pressure recovery as produced by a 
theoretical sharp lip loss mechanism. Reference 10 presents the theoretical derivations of 
Equation 144. The Mach number at the cowl lip, M u can be determined from continuity. 
Equation 145 is solved iteratively to find the inlet lip Mach number as a function of the inlet 
throat Mach number and the contraction area ratio from the inlet lip to the inlet throat. For 
cowl lips which are not sharp, but have some degree of bluntness, Equation 146 has been 
developed by the author to account for the effects of a non-zero cowl lip radius on the lip 
loss recovery given by Equation 144. Data from reference 11 was used to determine the 
exponential damping constant used in Equation 146. 


Results 


Results from the IP AC methodology are presented for three sample cases: a Mach 2.0 pitot 
inlet, a Mach 2.4 axisymmetric inlet, and a Mach 5.0 two-dimensional inlet. Example case 
output files, each containing copies of the respective input sets, can be found in Appendices 
II through IV. Additionally, a program User’s Guide which describes the input set and 
program usage can be found in Appendix I. 

The geometry of the pitot inlet sample case is shown in Figure 21. The pitot inlet is 
axisymmetric for this particular design and has a blunt cowl lip for improved low speed total 
pressure recoveries. Figure 22 shows a performance summary over the entire Mach number 
operating range for the inlet. The corrected airflow has been matched to a typical engine 
demand schedule, as shown in the top plot of Figure 22. The inlet throat Mach number was 
varied and used as an inlet control parameter in order to provide inlet/engine airflow 
matching. The resulting total pressure recovery and inlet drags are shown in the middle and 
lower plots in Figure 22. Note that the cowl drag is the dominant drag for this inlet design. 
This is an expected result of the blunt cowl lip feature of the inlet. 

Figure 23 shows the design and variable geometry features for the axisymmetric sample case. 
Both internal cowl surface variable geometry and a translating centerbody are used to control 
the operation of this inlet. Figure 24 shows a performance summary for the axisymmetric 
sample case over the entire Mach number range of inlet operation. Again, the inlet was 
designed and operated in accordance with a typical engine airflow demand schedule. The 
axisymmetric inlet is a mixed compression design with a starting Mach number of 1.6. This 
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inlet also requires a boundary layer bleed system. The sharp inlet lip can be seen to result in 
relatively lower take-off total pressure recoveries for this particular design. 

The inlet design and variable geometry features for the two-dimensional sample case is 
shown in Figure 25. This inlet uses a three ramp compression surface shock-on-lip design at 
Mach 5.0. The second and third ramps are movable and are used for inlet operation control. 
Figure 26 shows a performance summary for the two-dimensional design. The inlet employs 
both a boundary layer bleed system and an engine bypass system. The variable geometry 
ramp positions and bypass mass flow variations are used to provide matched airflow for a 
typical engine demand schedule. The inlet starting Mach number is 2.0 for this particular 
design. As is typical for high speed inlet systems, severe transonic drags are seen in the 
lower plot of Figure 26. 


Summary 

A series of analyses have been developed which permit the calculation of the performance of 
co mm on inlet designs. The methods presented are useful for determining the inlet weight 
flows, total pressure recovery, and aerodynamic drag coefficients for given inlet geometric 
designs. Limited geometric input data is required to use this inlet performance prediction 
methodology. The analyses presented here may also be used to perform inlet preliminary 
design studies. The calculated inlet performance parameters may be used in subsequent 
engine cycle analyses or installed engine performance calculations for existing uninstalled 
engine data. 
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List of Symbols 


A 

C D 

0)2 

Cp 

Ct 

D 

f 

Id 

8 

K 

K d 

K f 

K m 

K 0 

K a 

K, 

L 

L d 

M 

N 

P 
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r,R 

R 

T 

T t 

u 

v 

V 

W 

x , X 

y, y 

Yo 

cc 

& 

y 

Vv 

e 

0 

8 

P 


area, or cross-sectional flow area 
inlet capture area 
drag coefficient 

empirical cowl lip suction factor 
pressure coefficient 
oblique exit nozzle drag factor 
drag force, or diameter 

singularity distribution for the perturbation velocity potential 
diffuser friction factor 
gravitational constant 

empirical normal shock wave standoff factor 

empirical subsonic diffuser total pressure loss factor 

empirical subsonic diffuser friction loss factor 

empirical subsonic diffuser throat Mach number factor 

empirical subsonic diffuser offset loss factor 

empirical cowl lip suction factor 

empirical cowl lip suction factor 

inlet local location, or axial length 

diffuser axial length 

Mach number 

number of surface segments 

pressure 

total pressure 

dynamic pressure 

radial coordinate 

gas constant 

temperature 

total temperature 

axial perturbation velocity 

radial perturbation velocity 

velocity 

weight flow rate 

axial coordinate 

normal coordinate 

diffuser offset normal length 

forebody model angle 

shock wave angle 

ratio of specific heats 

discharge nozzle velocity coefficient 

referenced total temperature, or a flow/surface angle 

diffuser half-angle 

reference temperature for real gas model 

referenced total pressure 

density 
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X Mach number parameter in slender body theory 

v Prandtl-Meyer function 

£ integration parameter in slender body theory 

a correction factor for effective cowl angle 

<t> perturbation velocity potential, or a polar angle 

0 cowl curvature function 


Subscripts 

0 

1 

2 

ADD 

BL 

BLD 

BP 

BYP 

c 

C 

CWL 

cr 

e 

eff 

ENG 

i 

I 

L 

UP 

LS 

n 

N 

NS 

r 

S 

si 

SPL 

sub 

sup 

SY 

TH 

WAV 

x 

X 

* 

<t> 


free stream 
upstream, or cowl lip 
downstream, or engine face 
additive 

boundary layer bleed 
bleed 

engine bypass 
bypass 
cone 
cowl lip 
cowl 

critical operation 
effective 
effective 
engine flow 

general index, or ideal condition 

inlet capture 

inlet local 

inlet lip 

lip suction 

general index, or inlet ramp number 

normal component 

normal shock wave 

radial component 

centerbody or ramp surface 

sharp lip 

spillage 

sub-critical operation, or subsonic throat 
supersonic throat 

centerbody or ramp surface behind unstarted inlet normal shock wave 

inlet throat, or discharge nozzle throat 

wave 

axial component 
discharge nozzle exit 
sonic conditions 
polar component 
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Superscripts 


thermally perfect (non-ideal) gas property, or normalized velocity 
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Basic Inlet Modeling Elements 


u 

'c 

o 

U 1 

-Q 

D 

00 


c 

o 

*C0 

D 


C 

o 

'</) 

CO 

o £ 

c CL 

a> E 
~ o 
.E O 


c 

o 

’cn 

to 

O CD 

E fc. 
5 E 

X O 
LU CJ 


v co 
y O 
-c « 

CD «*— 
> LJ 




CD 

• - 0) 
cn o 
c a 
LJ L_ 



o 

o 



5 

o CL 

a 3 



^ o 
cd a 




E 

o 

<d a> 

<D l . 

it 00 


Figure 2 


Basic Inlet Modeling Elements Incorporated in the Analyses 
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Inlet Flow Stations 
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Figure 12 


Bleed System Operating Characteristics 
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Figure 13 Oblique Exit Nozzle Drag Factor 
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Figure 14 Effective Nozzle Discharge Pressure for Oblique Exits 
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Figure 15 Bleed and Bypass System Nozzle Exit Area Ratio 
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Figure 17 Cowl Lip Suction Factor, K a 
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Figure 18 Cowl Lip Angle Correction Factor, a 
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Figure 19 Cowl Lip Suction Factor, K & 
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Figure 20 Cowl Lip Suction Factor, C m 
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Figure 22 Pitot Inlet Performance 
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Axisymmetric Inlet Variable Geometry 
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IPAC - Inlet Performance Analysis Code 


Input List Description 


All variables are defined as implicit real*4 (a-h,o-z) unless 
otherwise noted in the following description. Variables 
beginning with the letters i-n are defined as integer unless 
otherwise noted. Any array variables are noted below with 
dimensions, ie . var(10). Default values are listed in the given 
assignments below. 


&ipac - namelist input set identifier, required 

table=' ipac . dat ' - tabular output data file name, character*80 


title=' ' 
echo=0 


iout=4*l 


f igure=0 


npts=10 , 20 


xmach0=0 . 01 
alt=0.. 0 
alpha0=0 . 0 
gama=l .4 


igas=0 


forbdy=0 


input case title, character*80 

echo flag, echoes input set to output if =1, 
integer 

output control flag array, setting each element of 
iout =1 writes additional data to output file, 
iout(l) program execution status messages 
iout (2) formatted performance summary pages 
iout (3) inlet flow station properties table 
iout (4) inlet geometry data summary 
iout (5) 

figure output flag, writes inlet figure data and 
output files if =1, integer 

number of points defining the engine face spinner 
or blunt cowl lip, npts(l), and subsonic diffuser 
contours, npts(2), when output is written using 
the figure=l option, npts(2) 

flight free stream Mach number 

flight altitude (ft) 

flight vehicle angle of attack (degrees) 

ratio of specific heats for atmosphere 

real gas effects flag, real gas calculations are 
performed if =1, typically only needed if xmachO 
is greater than 2.0 

vehicle forebody effects flag, no forebody if =0, 
initial conic forebody if =1, initial ramp 
forebody if =2, if =-l then forebody effects are 
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alphai=0 . 0 

xmlmO=l . 0 
ptlptO=l . 0 


idim=l 


ac=l . 0 


ar=l . 0 
ramps =0 

theta=0 . 0 


rleng=0 . 0 


xleng=0 . 0 


xcowl=0 . 0 
ycowl=l . 0 
cowls =0 

cowlth=0 . 0 


directly input through variables xmlmO and ptlptO, 
integer 

- array of forebody relative angles (degrees) used 
if forbdy =1 or 2, 

alphai(lO) [=lst_angle, 2nd_angle, . . .] 

- ratio of inlet local to free stream Mach numbers, 
used only if forbdy=-l 

- total pressure recovery ahead of inlet, used only 
if forbdy=-l 

- inlet type flag, symmetric 2-D pitot if =-l, 
axisymmetric pitot if =0, 2-D pitot if =1, 

2 -dimensional if =2, axisymmetric if =3, 
bifurcated 2 -dimensional if =4 

- inlet capture area (ft**2) , area will be computed 
if =-i and engine corrected weight flow data is 
supplied 

- inlet capture area aspect ratio, square or 
circular if =1 

- number of external 2-D inlet ramps (max 10) , or 
for an axisymmetric inlet conic centerbody set =1, 
integer 

- array of relative angles (degrees) of 2-D inlet 
ramps, or for an axisymmetric inlet conic 
centerbody set equal to the cone half -angle, 
theta (10) (=lst_angle, 2nd_angle, . . .] 

- array of radial lengths (ft) of 2-D inlet ramps, 
or the axisymmetric inlet conic centerbody length, 
do not use if the variable xleng is used, 
rleng(10) [=lst_length, 2nd_length, . . .] 

- array of axial lengths (ft) of 2-D inlet ramps, or 
the axisymmetric inlet conic centerbody length, do 
not use if the variable rleng is used, 

xleng (10) [=lst_length, 2nd_length, . . .] 

- cowl lip axial distance from inlet origin (ft) 

- cowl lip normal distance from inlet origin (ft) 

- number of segments defining the external cowl 
surface (max 10) , integer 

- relative angles (degrees) of external cowl 
surfaces, cowlth(10) [=lst_angle, 2nd_angle, . . .] 
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cowlrl = 0 . 0 


cowlxl 

rclip= 

a2ac=l 

xldd2= 

clof f = 

hubtip- 

thetac: 

nishck= 

xlipth= 

athac= 

xmth=l 

xmns=l . 
xtrans= 
aOac=l . 

athal=- 

ptrec=- 

ptrob=- 

ptreb=- 

ptrib=- 


- normalized radial lengths of external cowl 

surfaces, do not use if the input variable cowlxl 
is used, cowlrl(lO) [=lst_length, 2nd_length, . . . ] 

0.0 - normalized axial lengths of external cowl 

surfaces, do not use if the input variable cowlrl 
is used, cowlxl(lO) [=lst_length, 2nd_length, . . . ] 

.0 - normalized cowl lip radius, sharp lip =0 

. 0 - engine face flow area to inlet capture area ratio 

$.0 - subsonic diffuser axial length to engine face 

diameter ratio 

).0 - normalized inlet origin to engine face centerline 

offset distance 

=0.3 - engine face spinner to fan tip radius ratio 

=0.0 - cowl lip internal angle (degrees) 

=-l - number of inlet internal shock wave reflections, 

calculate number of shocks if =-l 

=-l - normalized length of inlet internal duct from cowl 

lip to throat, calculate length if =-i 

1 - inlet throat area to inlet capture area ratio, 

calculate if =-l 

3 - inlet throat Mach number, calculate throat Mach 

number if =-l 

35 - inlet supercritical normal shock Mach number 

=0.0 - normalized centerbody translation distance 

0 - stream tube capture area ratio, usually calculated 

and not used as an input variable 

1 - inlet throat area to cowl lip flow area ratio, 

calculate if =-l 

1 - inlet total pressure recovery, calculate if =-l 

1 - total pressure recovery across oblique shock 

waves, calculate if =-l 

1 - total pressure recovery across external oblique 

shock waves, calculate if =-i 

1 - total pressure recovery across internal oblique 

shock waves, calculate if =-l 
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ptrns=-l 


ptrfr=-l 


ptrdf =-l 


ptrlp=-l 


fd=0 .0025 
bleed=0 . 0 


pblpt0 = 0 . 0 


thexbl=15 . 0 


nvbl=0 . 98 


nozzbl=l 


axthbl=l . 0 


bypass=0 . 0 


pbppt2=0 . 0 


thexbp=15 . 0 


total pressure recovery across normal shock wave, 
calculate if =-l 

total pressure recovery factor resulting from 
inlet surface friction ahead of the throat, 
calculate if =-l 

total pressure recovery factor resulting from 
subsonic diffuser behind the inlet throat, 
calculate if =-l 

total pressure recovery factor resulting from cowl 
lip flow losses, calculate if =-l 

subsonic diffuser friction loss factor 

array of inlet bleed flow mass fractions for each 
bleed system, up to 10 separate bleed systems can 
be defined, 

bleed (10) [=lst_sys_f rac, 2nd_sys_f rac , . . .] 

array of total pressure recovery in bleed plenum 
to freestream for each separate bleed system, 
pblptO (10) [=lst_sys_rec, 2nd_sys_rec, . . .] 

array of bleed flow discharge angles (degrees) 
relative to freestream for each bleed system, 
thexbl (10) [=lst_sys_angle , 2nd_sys_angle , . . . ] 

array of bleed flow discharge nozzle velocity 
coefficients for each separate bleed system, 
real*4 nvbl(10) [=lst_sys_coef , 2nd_sys_coef , . . . ] 

array of the type of bleed flow discharge nozzle 
used for each separate bleed system, convergent 
nozzle if =1, convergent -divergent nozzle if =2, 
nozzbl (10) [=lst_sys_type,2nd_sys_type, . . .] 

array of bleed flow discharge nozzle exit area to 
nozzle throat area ratio for each separate bleed 
system, set =1 if the nozzle is convergent, 
axthbl (10) [=lst_sys_ratio, 2nd_sys_ratio, . . .] 

array of inlet bypass flow mass fractions for each 
bypass system, up to 10 separate bypass systems 
can be defined, 

bypass (10) [=lst_sys_f rac, 2nd_sys_frac, . . .] 

array of total pressure recovery in bypass plenum 
to engine face for each separate bypass system, 
pbppt2 (10) [=lst_sys_rec, 2nd_sys_rec, . . .] 

array of bypass flow discharge angles (degrees) 
relative to freestream for each bypass system. 
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nvbp=0 . 98 


nozzbp=l 

axthbp=l . 0 

cdcowl=-l 

refcd=-l 

etype=0 

escale=l . 0 
fn=0 . 0 

sfc=0 . 0 

w2cor=0 . 0 

w2abs=0 . 0 

pt8pt2=l . 0 

refrec=-l 


thexbp(lO) [=lst_sys_angle, 2nd_sys_angle , . . .] 

- array of bypass flow discharge nozzle velocity 
coefficients for each separate bypass system, 
real*4 nvbp(10) [=lst_sys_coef , 2nd_sys_coef , . . . ] 

- array of the type of bypass flow discharge nozzle 
used for each separate bypass system, convergent 
nozzle if =1, convergent -divergent nozzle if =2, 
nozzbp (10) [=lst_sys_type, 2nd_sys_type, . . .] 

- array of bypass flow discharge nozzle exit area to 
nozzle throat area ratio for each separate bypass 
system, set =1 if the nozzle is convergent, 
axthbp(lO) [=lst_sys_ratio, 2nd_sys_ratio, . . . ] 

- cowl drag coefficient, sum of lip and wave drags, 
calculate if =-l 

- reference inlet drag coefficient, will be set 
equal to -cdcowl if =-l 

- array of engine type for each engine in an engine 
module, up to 10 engines per module, set =1 for a 
ramjet engine, set =2 for a turbojet engine, 
integer etype(lO) [=lst_eng_typ, 2nd_eng_typ, . . . ] 

- array of sizing scale factors for each engine, 
escale(lO) [=lst_eng_size, 2nd_eng_size, . . .] 

- array of the uninstalled net thrust (lb) for each 
engine in an engine module, 

fn(10) [=lst_eng_thrust , 2nd_eng_thrust , . . .] 

- array of the uninstalled specific fuel consumption 
(lbm/hr/lbf) for each engine in an engine module, 
sfc(10) [=lst_eng_sfc, 2nd_eng_sfc, . . .] 

- array of the uninstalled engine face corrected 
weight flow (lb/s) for each engine, 
w2cor(10) [=lst_eng_flow, 2nd_eng_f low, . . .] 

- array of the uninstalled engine face absolute 
weight flow (lb/s) for each engine, 

w2abs (10) t=lst_eng_f low, 2nd_eng_f low, . . . ] 

- array of the total pressure ratio across the 
engine, from nozzle throat to engine face, 
pt8pt2 (10) [=lst_eng_ratio, 2nd_eng_ratio, . . . ] 

- array of the reference total pressure recovery 
used for each engine, set = -l for MIL-SPEC, 
refrec ( 10 ) [=lst_eng_rec, 2nd_eng_rec, . . . ] 
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nozzle 


noeng=l 

aero 


&end 


array of engine module nozzle data, real*4 
nozzle (1) uninstalled engine data Cfg 
nozzle (2) actual nozzle gross thrust coefficient 
nozzle (3) actual nozzle drag coefficient 
nozzle (4) reference area (ft**2) for nozzle Cd 


number of engine modules on vehicle 


array of 
aero (1) 
aero (2) 
aero (3 ) 
aero (4) 


the flight vehicle aerodynamic data 
lift coefficient 
drag coefficient 
angle of attack 

reference area (ft**2) for Cl and Cd 


namelist identifier, required 


Notes on Input Usage 

The input and output filenames may be specified on the command 
line after the program name. The extensions .in and .out may be 
left off the filenames and will automatically be appended. 

system_prompt> ipac ipac.in ipac.out 

The program IPAC reads the namelist input set from an input file 
(the default is ipac.in) and executes the required calculations 
for that case. The output is written to an output file (the 
default is ipac.out) and to another tabular data file specified 
by the input variable table in the namelist input set. If there 
are subsequent namelist input sets in the file, they in turn are 
executed, and in this manner numerous cases can be run to design 
and/or analyze an inlet system over a range of operating 
conditions. Since the program uses namelist input reads, if a 
variable is defined once in an input set, it is not necessary to 
redefine it again in subsequent input sets, unless the value 
changes. Also, since nearly all of the input variables have 
predefined defaults, it is usually only necessary to assign 
values to a few variables to run the program properly. The 
character string pairs /* ...comments... */ are parsed and 
discarded by the input set read routine, thus allowing for the 
inclusion of comments, or the exclusion of commented out inputs, 
in the input file. 

There are a few subtleties which the user needs to be aware of to 
effectively use IPAC. The following paragraphs describe some of 
the ways the various input variables are used to model inlet 
systems . 

General Output Control : The first 6 variables listed above 

determine the output features for IPAC. The data file defined by 
the table variable will contain a summary tabular dataset of 
inlet operation and performance quantities such as: pressure 
recovery, mass flow ratios, and drag coefficients. These 
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quantities are sufficient to compose a set of inlet performance 
maps. To facilitate the generation of performance maps, more 
than one data file can be defined by the table variable in 
subsequent namelist input sets. Thus, a range of inlet operating 
points can be written to different tabular datasets. The user 
must then re -format these datasets to construct inlet map files 
appropriate for other analysis codes. 

The title variable is printed for each output case if defined. 

The echo variable can, and is recommended, to be set to 1. This 
will print the namelist input set ahead of each output case. 
Additionally, if echo is set to 2 then the entire input file will 
be printed at the top of the output file. The array variable 
iout is used to control the level of data written to the output 
file, ipac.out. Setting the elements of iout =1 will result in 
additional output data. Currently there are 4 elements in iout 
which can be used for output control . Status messages of program 
execution information are enabled/disabled by iout(l) =1/0. 

These single line printouts of pertinent variable values from 
each major analysis segment (as the code executes) are useful for 
quickly assessing the progress of the inlet design, operation, 
and performance modeling. Printout of formatted inlet 
performance summary data is enabled/disabled by iout (2) =1/0. A 
formatted data table of flow properties at each of the inlet flow 
stations is enabled/disabled by iout (3) = 1/0. A brief inlet 
geometry data summary is enabled/disabled by iout (4) = 1/0. The 
program defaults will print all of the above information for each 
input case . Complete inlet performance data is written to 4 
other tabular datasets *.dat for all input cases executed. This 
information is very easily graphed by a plotting package of the 
user's choice. 

Printout of the inlet geometry contours is enabled by setting the 
figure variable =1. Additional output files *.fig are written 
which contain (x,y) coordinate pairs that can be used to 
construct a simple line drawing of the inlet geometry, and which 
can be viewed by the user's own plotting package of choice. 

The figure variable should be set to 1 in only one input set, and 
then reset to 0 for the rest of the cases since the *.fig output 
files are overwritten for each case. The array input variable 
npts can be used to increase the number of points written which 
define. the subsonic diffuser, blunt cowl lips, and engine face 
segments of the figure. This allows for greater resolution of 
the curved surfaces in the geometry. 

Flight Conditions : The Mach number and altitude for flight are 

set in variables machO and alt. If a positive number is 
assigned to alt then the program will use that value for the 
altitude in ft. If alt is assigned a negative number, then the 
program will assume that the user has entered a flight dynamic 
pressure (in psf) instead, and will find an appropriate altitude 
for the specified flight Mach number. This is a convenient 
feature for finding constant Q flight paths. If the vehicle is 
situated at an angle of attack to the freestream, the variable 
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alphaO should be used. If the user feels it is necessary to 
adjust the ratio of specific heats constant for the atmosphere, 
the variable garna can be used. If flight conditions exceed Mach 
2.0, it is recommended that igas be set to 1 to adjust ideal gas 
assumptions for real gas effects which become important for high 
speed flight . 

Vehicle Effects : If the inlet is located close to the body /wing 

of the vehicle it may be necessary to account for changes in flow 
conditions entering the inlet as a result of vehicle effects. 

The variable forbdy controls how the vehicle effects are modeled. 
Values of 1 or 2 assigned to forbdy can model simple combinations 
of conic and ramp configurations. The necessary relative angles 
(degrees) are input through the array variable alphai. 

Compressive turning is denoted by a positive angle, and 
expansions are denoted by a negative angle. If a very complex 
flowfield is produced by the vehicle, the changes in Mach number 
and total pressure can be directly input in variables xmlmO and 
ptlptO (provided these values are known) if forbdy is set to -1. 

Inlet Geometry : A number of variables are used to describe the 

inlet geometry to be modeled. The first is idim which specifies 
the basic inlet type: pitot, axisymmetric, or 2-dimensional. The 
permitted values of idim follow. 


idim =- 1 
0 
1 
2 

3 

4 


symmetric 2-D pitot inlet 
axisymmetric pitot inlet 
2-D pitot inlet 
2 -dimensional inlet 
axisymmetric inlet 
bifurcated 2 -dimensional inlet 


If the inlet is 2-D the aspect ratio, variable ar, is the inlet 
width divided by height. If the inlet is axisymmetric then ar is 
interpreted as fraction of a full-circle. Thus, for a hemi- 
circular axisymmetric inlet, ar would be set to 0.5. The 
variable determining the gross size of the inlet is the capture 
area, ac in square ft. This can be simply set to 1 for easy 
normalizations, any physical size in square feet, or if set to -1 
will be calculated and automatically sized to match the engine 
demand airflow requirements if this data is supplied. 

External Compression Surfaces: The variables ramps, theta, 

rleng, and xleng define the inlet external compression surfaces 
for axisymmetric and 2-D inlets. For axisymmetric inlets, ramps 
must be set to 1, and theta is set to the conic centerbody half- 
angle. Either rleng or xleng, in ft, can be used to define the 
centerbody length, but not both. For 2-D inlets, ramps can be 
set up to a maximum of 10, and theta is then set to the relative 
angles (degrees) of each ramp. Either rleng or xleng can be used 
to define the lengths of each ramp, but not both. It is 
recommended that rleng be used since it does not change as the 
ramp angles are varied. 
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Cowl Lip & Shock-On-Lip Design Feature: The location of the cowl 

lip is specified by variables xcowl and ycowl in ft. These 
variables are used in both axisymmetric and 2-D inlets. For 
axisymmetric inlets ycowl is the radial distance from the inlet 
centerline. There is a feature in IPAC which will automatically 
calculate the location of the cowl lip for the shock-on-lip 
condition. Also, this feature will calculate the ramp lengths 
for multiple ramp 2-D inlets, placing all of the shock waves on 
the cowl lip, provided that the ramp angles are specified. This 
is a very useful design feature. To use this automatic design 
capability do the following in the very first namelist input set. 

(1) set ramps to the number of ramps or 1 for a centerbody 

(2) set theta to the ramp or centerbody relative angle (s) 

(3) set rleng and xleng to 0.0, this is the program default 

(4) set xcowl =0 and ycowl =1, also the program default 

IPAC will then calculate the location of the cowl lip, and the 
lengths of all the ramps for shocks-on-lip for the specified 
flight Mach number, xmachO . These results will be remembered for 
subsequent cases, and there is no need to input these values by 
hand. 

External Cowl Surfaces: The variables cowls, cowlth, cowlrl, and 

cowlxl define the external contour of the inlet cowl surface. 

The number of segments is specified in cowls, the relative angles 
(degrees) in cowlth, and the lengths in either cowlrl or cowlxl. 
The lengths are normalized by ycowl and thus specified as 
multiples of ycowl. A blunt cowl lip radius can be specified by 
the variable rclip and this radius is also normalized by the 
length ycowl . 

Subsonic Diffuser: There are 4 input variables which are used to 

define the geometry of the subsonic diffuser element in an inlet. 
The engine face flow area is defined as a ratio relative to the 
inlet capture area through the variable a2ac. The axial length 
of the diffuser is defined as a ratio relative to the engine face 
diameter through the variable xldd2 . The vertical offset 
location of the engine face is defined as a normalized distance 
from the inlet origin to the engine centerline, through the 
variable cloff, as a multiple of the distance ycowl. The 
variable hubtip performs a number of functions. If hubtip is a 
positive number then it defines the engine face spinner to fan 
tip radius ratio. If hubtip equals 0.0 then no engine spinner 
exists but the engine face is still assumed to be circular. If 
hubtip is a negative number then the program will recognize that 
the user has indicated that the engine face is not circular, but 
rather 2 -dimensional , and that the value specified in hubtip is 
now the aspect ratio for the 2-D engine face duct area. 

Internal Shocks : For supercritical operation of mixed 

compression inlets, internal shock waves are formed between the 
cowl lip and the inlet throat. The model used in IPAC is 
relatively simple for this internal supersonic duct. A 
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constantly converging channel is used to model the flow from 
inlet cowl lip to throat regions. The difference between the 
internal cowl lip angle (degrees) , thetac, and the last external 
ramp angle forms the net convergence angle for the duct model. A 
single shock wave train, reflecting off each duct wall, is used 
to model the supersonic flow. The variable nishck is used to 
specify how many shock waves will be permitted in the duct, and 
this value will be calculated if set to -1. The variable xlipth 
is the normalized length (multiple of ycowl) of the duct from the 
cowl lip to the throat, and will also be calculated if set to -1. 
The variable athac is the inlet throat area to capture area 
ratio. This variable is critical in determining the inlet 
operation. If athac is set to -1 this ratio will be calculated. 
The variable xmth is the inlet throat Mach number. By specifying 
an inlet throat Mach number and area, the mass flow of the inlet 
is uniquely determined. 

In a typical design point calculation it is easiest to specify 
the throat Mach number, xmth, and then for supercritical 
operation the rest of the variables, nishck, xlipth, and athac 
will be determined. For subsequent calculations, the inlet 
throat area will then be determined from the inlet geometry, and 
the throat Mach number will in turn be calculated. The variable 
■ vmn a is the Mach number ahead of the internal terminal normal 
shock. Note that for supercritical operation xmns must be 
greater than the throat Mach number xmth. As the normal shock 
Mach number is increased, the shock will be pulled further 
downstream from the inlet throat into the subsonic diffuser. 

This will also decrease the inlet recovery. Specifying xmns is 
another control variable which can be used to match the inlet 
supply corrected airflow to the engine demand. If xmns is set =0 
and the inlet is operating supercritical, then the flow at the 
engine face will be calculated as supersonic flow. This permits 
the modeling of supersonic through-flow fan and scramjet inlets. 

Variable Geometry : After the inlet design point is calculated in 

the first namelist input set, the throat area can be increased or 
decreased by variable geometry features for off -design operation. 
For multi-ramp inlets, the ramp angles theta can be redefined by 
the user in subsequent namelist input sets. The cowl internal 
angle thetac can also be changed. A very common variable 
geometry mechanism for axisymmetric inlets is the translating 
centerbody, and the input variable xtrans can be used to move the 
centerbody forward a specified distance which is a multiple of 
ycowl. Note that xtrans works only for axisymmetric inlets, and 
produces no translation for two-dimensional inlets. 

Two additional, although not typically used, input variables are 
the stream tube capture area ratio, aOac, and the throat to cowl 
lip flow area ratio, athal. The stream tube capture area ratio 
is usually calculated by the program, however, it is possible 
that for some inlets the capture area ratio can be defined, and 
then for a given geometry the inlet throat Mach number would be 
calculated . 
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Recovery Overrides : All of the input variables beginning with 

ptr are the total pressure ratios for various loss producing 

mechanisms and are normally calculated in the program. The user 
has the option of overriding these calculations and directly 
entering values for any and all of these terms. Normally this is 
not done, however, if other more complex analyses have been 
performed for an inlet design, then the user can use those values 
instead of the ones that IPAC would normally calculate. 

An additional input variable is the friction loss factor, fd, 
which is used in the subsonic diffuser loss model. The default 
value is 0.0025 and this value is suitable for most typical 
subsonic diffuser designs. 

Bleed and Bypass Systems : Boundary layer bleed is a necessary 

component for all high speed inlet systems. In order to 
stabilize the shock wave boundary layer interactions, a small 
amount of air is removed through the walls of the inlet. This 
air is then dumped overboard and a momentum drag is incurred. 

Mass removed and dumped ahead of the inlet throat is called 
bleed, and is necessary for inlet operation. Mass removed and 
dumped behind the throat is called bypass, and is sometimes 
necessary for inlet/engine matching. Up to ten independent bleed 
and ten independent bypass systems can be defined. Both bleed 
and bypass inputs work the same way, and therefore, only the 
bleed variables will be directly discussed. The user must 
specify the variable bleed, the fraction of captured airflow 
which is to be dumped. The variable pblptO is the total pressure 
ratio (bleed plenum to freestream) for the bleed system and must 
also be chosen. The rest of the variables, thexbl, nvbl, nozzbl, 
and axthbl may be left at the default values. For bypass 
systems, since there is typically much more pressure available 
for expansion, a convergent -divergent nozzle may be used. 

The input variable bleed can be defaulted to any negative number 
to automatically calculate the amount of boundary layer bleed as 
a function of inlet local Mach number. If bleed set =-l, then a 
single bleed system will use the default bleed rate. If bleed is 
set =-0.8, then a single bleed system will use 80% of the default 
bleed rate. If bleed is =-1.5, then a single bleed system will 
use 150% of the default bleed rate. If bleed =-0.4, -0.5, then 
two bleed systems will use a total of 90% of the default bleed 
rate. The bleed plenum total pressure recovery variable, pblptO, 
can also be defaulted to an internal calculation, again as a 
function of inlet local Mach number. Set pblptO to: -1 for the 
nominal average recovery, -2 for the high pressure porous 
recovery, -3 for the low pressure porous recovery, and -4 for the 
throat slot recovery. Each individual bleed system can use any 
appropriate bleed configuration recovery. 

For the bypass system, the total pressure recovery in the bypass 
duct, pbppt2, can be calculated from a bypass duct loss as a 
function of bypass fraction. To calculate, set pbppt2 =-l. 
Typically, when matching inlet supply and engine demand the inlet 
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provides excess airflow which must be bypassed. If engine data 
is supplied, and the inlet has excess airflow capacity, set 
bypass =-l to automatically match the inlet and engine airflows 
by adjusting the bypass fraction. This feature only works on the 
first bypass system, the other bypass systems if defined cannot 
be automatically matched but must be directly input. 

Drag Accounting : The exact details of which inlet drag 

components should be charged to propulsion or airframe are a 
subject of continual debate. Most notable is the cowl drag, 
which is comprised of cowl blunt lip and cowl wave drags. IPAC 
calculates these drag components if the input variable cdcowl is 
set = -l, the program default. If any other positive value is 
assigned to cdcowl that value will be used, and the lip and wave 
drag calculations will be skipped. Since the external drag on an 
engine nacelle is often accounted for in the vehicle aerodynamic 
performance, another input variable refcd has been included. 

This variable represents the reference drag coefficient for the 
inlet installation, and thus part of the inlet drag can be 
accounted for in the vehicle aerodynamic data. 

The net inlet drag at an engine operating point is called the 
power setting drag, and the power setting drag is equal to the 
total of all the inlet drags (spillage, bleed, bypass, cowl lip 
and wave) less the reference drag. Often the cowl drag 
components are accounted for in the vehicle aerodynamic data. If 
refcd is set =-l, then the reference drag will be set equal to 
the inlet cowl drag. This will result in an inlet power setting 
drag comprised of only spillage, bleed, and bypass drags. This 
is the default for the program, where refcd is =-l. 

Engine Data : Engine data can be supplied to the program, and 

IPAC will perform installation calculations and re-calculate 
engine data if desired. It is assumed that engines are in 
separate modules, and there can be more than one engine (up to 
10) in a module. However, each module has a single inlet, and 
possibly a common nozzle. The variable etype specifies the types 
of engines in a module. The only types are turbojet and ramjet 
at this time. A turbofan can be modeled as two separate 
turbojets, one with and one without fuel. The variable escale 
can be used to adjust the size of the engines for inlet/engine 
matching and sizing studies. 

Each engine in the engine module is specified as an element in 
the array variables: fn, sfc, w2cor, w2abs, pt8pt2, and refrec. 
The engine net thrust and specific fuel consumption are specified 
by the input variables fn and sfc. To perform the installation 
calculations the total pressure ratio across the engine, pt8pt2, 
and the inlet recovery used in determining the uninstalled engine 
data, refrec, must also be supplied. The user may specify that a 
MIL-SPEC inlet recovery was used in the., uninstalled engine data 
by setting refrec =-l, which is the program default. The 
absolute engine weight flow must be specified in the input 
variable w2abs in lb/s. The corrected engine airflow is also 
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required in the variable w2cor in lb/s. To perform engine 
installation calculations correctly, the uninstalled engine 
corrected weight flow must be the same as the inlet supply 
corrected weight flow. 

The equality of inlet supply and engine demand corrected weight 
flow is called inlet /engine matching. All proper inlet designs 
must be matched to an engine demand corrected airflow schedule. 

If the user is designing an inlet and there is no engine data 
available, the program will construct an engine demand corrected 
weight flow schedule automatically for a "typical" engine. If 
w2cor is set =-l at the inlet design point, then the program will 
automatically calculate an engine demand corrected weight flow 
which matches the inlet supply corrected weight flow. By leaving 
w2cor =-l for the rest of the inlet operating points, the program 
will calculate the "typical" engine demand corrected weight flow 
schedule as a function of flight Mach number. The user may then 
use this schedule of engine demand corrected weight flow for 
inlet/engine matching over the off-design operating points. 

The installed thrust for the engine module will be equal to the 
uninstalled engine thrust adjusted for the actual inlet recovery 
less the inlet power setting drag. Note that engine data that is 
installed but not properly matched with the inlet supply 
corrected weight flow is fundamentally incorrect since 
conservation of mass will be violated. 

Nozzle Data : If nozzle data is available, the installation 

calculations will also adjust the engine data for nozzle effects. 
The inputs are in the array variable nozzle, and include the 
gross thrust coefficient used in the engine data, the actual 
gross thrust coefficient for the nozzle used, a drag coefficient 
for the nozzle, and a reference area. Note that when using the 
nozzle input variable it is assumed that only one nozzle is used 
for each engine module, even though more than one engine can be 
in a module . 

Vehicle Data : Since it is often of interest to see how engine 

systems size on the vehicle, IPAC can accept vehicle aerodynamic 
data. Thus engine sizing studies can also be performed. The 
variable noeng sets the number of engine modules on the vehicle. 
The array variable aero contains the vehicle lift and drag 
coefficients, angle of attack, and reference area. Thus, the 
program can install engines with an inlet design, and can then 
determine if the propulsion system is capable of powering the 
aircraft throughout the flight regime. 
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Mach 2.0 Pitot Inlet 
Example Case 
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Figure II. 2 


Mass Flow Ratios 
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Figure II. 4 
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Figure II.5 Corrected Airflows 
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Appendix IV 


Mach 5.0 Two-Dimensional Inlet 
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Figure IV. 1 Total Pressure Recoveries 
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Figure IV. 3 Drag Coefficients 
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Figure IV. 5 Corrected Airflows 
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Inlet Mass Flow* 
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Inlet Flow Properties free inlet cowl throat engine 

stream local lip face 
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